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ABSTRACT 

Most millisecond pulsars, like essentially all other radio pulsars, show timing errors well in excess 
of what is expected from additive radiometer noise alone. We show that changes in amplitude, 
shape and pulse phase for the millisecond pulsar J1713+0747 cause this excess error. These changes 
appear to be uncorrelatcd from one pulse period to the next. The resulting time of arrival variations 
are correlated across a wide frequency range and are observed with different backend processors on 
different days, confirming that they are intrinsic in origin and not an instrumental effect or caused 
by strongly frequency dependent interstellar scattering. Centroids of single pulses show an rms phase 
variation « 40 /zs, which dominates the timing error and is the same phase jitter phenomenon long 
known in slower spinning, canonical pulsars. We show that the amplitude modulations of single 
pulses are modestly correlated with their arrival time fluctuations. We also demonstrate that single- 
pulse variations are completely consistent with arrival time variations of pulse profiles obtained by 
integrating TV pulses such that the arrival time error decreases proportional to 1/ ^/N . We investigate 
methods for correcting times of arrival for these pulse shape changes, including multi-component 
TOA fitting and principal component analysis. These techniques are not found to improve the timing 
precision of the observations. We conclude that when pulse shape changes dominate timing errors, 
the timing precision of PSR J1713-I-0747 can be only improved by averaging over a larger number of 
pulses. 

Subject headings: gravitational waves, methods: statistical, pulsars: general, pulsars: specific (PSR 
J1713+0747) 
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1. INTRODUCTION 

The timing of pulsars has enabled many studies of funda- 
mental importance to physics and ast rophysics, including 
preci s ion tests of general relativity (jTavlor k, Weisbera 
119821 : Kramer et al.l I2006D . strong constraints on nu- 
clear equations of s tate (jLattimer fc PrakashI 120071 : 
[Demorcst ct a]] l2010f ). and the discovery of the 
first Eart h-mass planets outside of the solar system 
(fWols zczan fc Fraill [19921 ) . By studying the correlated 
variations in pulse times of arrival (TOAs) from a set 
of millisecond pulsars (MSPs) that comprise a pulsar 
timing array (PTA), it is poss ible to detect nanohertz 
gravitational waves (GWs) ( HelUngs fc Downs! [l98l 
[Foster fc Backedll990[ ). Candidate sources of GWs in the 
nanohertz frequency band are binary supermassive black 
holes in the throes of merger, o scillating cosmic s trings, 
and the inflation era Universe penet et al.ll2009() . The 
net perturbation due to these sources to TOAs is ex- 
pected to be at most tens of nanoseconds over 5—10 
years. 

The most plausible source of gravitational waves is a 
stochasti c background associated with massive blac k hole 
binaries ()Jaffe fc Backed 120031: iSesana et al.ll2008[ ). For 
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this background, a modestly significant (< Aa) detec- 
tion can be ac hieved by observin g 20 to 40 pulsars over 
5 to 10 years (jJenet et al.l 120051 ) assuming monthly ob- 
servations with 100 ns timing precision. 

The required timing stability of < 100 ns can only be 
achieved with MSPs, pulsars with spin periods between 
« 1 and 20 ms, which are more spin stable than canonical 
pulsars (CPs, pulsars with spin periods between w 30 ms 
and 8 s). Over periods of years, CPs show large levels of 
correlated spin noise likely associated with instabilities 
in bot h the pulsar magnetosp here and th e neutron star 
itself (ICOTd es fc D o^ [19851 iD^A lcssandro et al.l [19951: 
iLvne et al.ll2010E ,Shannon fc Cordesi,2010 ^ 

The precision of pulsar timing is at minimum limited by 
radiometer noise associated with the telescope receiving 
system and sky background. If this is the only form 
of noise, the sensitivity of observations is improved by 
increasing observation bandwidth, increasing observing 
time, or by increasing telescope sensitivity through larger 
collecting area or lower system temperature. 

However, there are certainly other processes that affect 
the timing precision of pulsars and limit the sensitivity 
of PTAs to GWs. The sources and i mportance of a wide 
range of effects are summarized in iCordes fc ShannonI 
(|2010( ). Stochastic timing variations can be broadly bi- 
furcated into those with stationary statistics, such as 
white noise, and others that have or appear to have 
nonstationary statistics, such as those with very broad 
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fluctuation spectra, including steep power laws with 
most of the power at low frequencies, often described 
as "red" processes. Red-like timing noise appears to 
describe the residuals seen in many canonical pulsars 
though it is still under debate whether on the longest 
times scale, the timing variations have non-stationary 
stati stics or show band-limited or quasi-period ic statis- 
tics HLvne et al.ll2010l: [Shannon fc Cordesll20Tot) . Propa- 
gation of radio ernission through the interstellar medium 
(lArmstrondHOSl iColes et al.llMol : iCordes fc Shannon! 
120101 ) also produces timing variations that have red 
power spectra. 

On shorter time scales ranging from one pulse period to 
integrations from several minutes to an hour, TOA vari- 
ations exceed what is expecte d from radiomete r noise 
alone. In a study of CPs, ICordes fc Downs! (|1985f ) 
demonstrated that timing errors were generally always 
larger than expected from additive radiometer noise and 
interpreted the excess in terms of changes in pulse shape 
and amplitude on pulse-to-pulse time scales. Termed 
"jitter," this phenomenon has been seen in virtually all 
pulsars that have been closely investigated. Integrated 
pulse profiles include these variations and thus so do the 
TOAs derived from them. 

While the link between single pulses and timing vari- 
ations is well established in CPs, studies of MSPs are 
fewer because MSPs typically have much lower flux den- 
sities. However, a few cases indicate that single-pulse 
variations in MSP s are consistent with those seen in CPs. 
iJenet et all ()1998D studied the variability of single pulses 
of the bright millisecond pulsar J0437— 4715. They found 
that single pulses contained a wide variety of morpholo- 
gies and concluded the emission showed similar proper- 
ties to that observed in canonical pulsars. Oslowsk i et al.! 
(|201l!) followed up this analysis by identifying a correla- 
tion between arrival time and pulse shape changes, and 
used this correlation to improve timing precisio n by 20%. 
In a m ore recent study of PSR J0437— 4715, iLiu et all 
(|2012l ) found excess tim ing errors that they a t tribut ed to 
pulse shape variations. Edwards fc Stappers (2003) ana- 
lyzed the intensity modulation of a set of northern MSPs 
and suggested that there are periodic intensity modula- 
tions in MSP arrival times comparable to the drifting 
sub- pulse phenomenon observed in many canonical pul- 
sars. IJenet fc Gill (poM) found remarkable intensity sta- 
bility of the main pulse of PSR B1937+21, but accom- 
pani ed by large intensity m odulation in giant pulse emis- 
sion. iKramer et al.! (|1999l ) showed that the relative inten- 
sity of the components of the MSP J1022-I-1001 change 
on longer time scales and that a better timing solution for 
this pulsar could be obtained when these changes were 
incorporated. 

In this paper, we focus on high S/N observations of 
the millisecond pulsar J1713-I-0747 to analyze fast pulse 
shape variations through measurements of single pulses 
and average pulse profiles. In Section [21 we summarize 
the observations and the data reduction procedure. In 
Section [31 we infer both directly and statistically that 
intrinsic pulse shape changes cause TOA variations in 
this pulsar. In Section [H we rule out other processes 
that mimic pulse jitter in our observations. In Section [5l 
we present two different correction procedures that can 



be used to correct for TOA variations and apply the tech- 
niques (unsuccessfully) to our observations. In Section [6l 
we discuss observing strategies that mitigate the effects 
of pulse profile variations on timing precision. 

2. OBSERVATION AND ANALYSIS PROCEDURE 

The relatively nearby MSP J1713-h0747 has a period P = 
4.57 ms, dispersion measure D M= 16 pc cm~^, and a 
pulse width of 110 /is (FWHM) ( Manchester et al.ll2005l ). 
Its large (period averaged) flux density of 7.4 mJy at 
1.4 GHz contributes to its being one of the best MSPs in 
terms of timing precision. We observed PSR J1713-I-0747 
at 1.5 GHz with the Arecibo 305-m telescop43 using the 
"L-band Wide" receiver, which provides two channels, 
one for each hand of circular polarization. The results 
presented here use data recorded with two backend signal 
processing instrume nts: the Wide-ban d Arecibo Pulsar 
Processor (WAPPs, iDowd et al.l[2000[) and the Arecibo 
Signal Processor (ASP. lDemorestll2007l ) . The ASP back- 
end is used for several long-term timing programs, in- 
cluding precision timing for GW detection. We use the 
WAPP data to study single-pulses but the ASP data to 
analyze TOA and pulse-shape variations of average pro- 
files. The combination of data then allows us to connect 
single-pulse results to the long-term timing programs on 
PSR J1713-f 0747 that make use only of average profiles. 

The WAPP data provided ~ 50A: pulses over a 4-minute 
time span on MJD 54983 and were used to analyze single 
pulses and sub-integrations of up to 512 pulses. Specif- 
ically, the WAPP instrument outputs a time series of 
autocorrelation functions (ACFs) of both polarizations 
with 3-level sampling. Offline, the sigproc packagc0 was 
used to Fourier transform the ACFs into spectra with 
192 frequency channels across a 100 MHz total band- 
width spanning 1.57 to 1.67 GHz with 32 //s resolution. 
After excising RFI, the data were dedispersed to form 
a single time series in the sum of the two polarization 
channels, which we call t he total intensity. We used the 
dspsr software package (jvan Straten fc Baiie^lMl to 
fold the data into subintegrations containi ng 1 to 1000 
pulse s, and the psrchive software package (jHotan et al.! 
12004! ) to measure TOAs fr om the subintegrations. 

The ASP backend removes dispersion coherently to pro- 
duce dedispersed time series in all four Stokes parameters 
which were then synchronously averaged in real time us- 
ing the pulsar ephemeris. Averages were outputted every 
10 s for 16 sub-bands, each with 4 MHz bandwidth, for 
a total of 64 MHz centered at 1410 MHz. 

The ASPFitsRea der pipeline was used to calibrate and 
generate TOAs (jFerdman! 120081 ) from ASP measure- 
ments. Polarization calibration and absolute flux calibra- 
tion were completed by comparing a pulsed signal gen- 
erated by a noise diode while the telescope was pointed 
on and off a calibrating source to the pulsed signal while 
the telescope was pointed at the pulsar. The radio galaxy 
CTD 93 was used as the calibration source because it is 

The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which was operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 

* sigproc.sourceforge.net 
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known to be both unp olarized and to have no observed 
flux variability ('Shaff er et al J 119991 ) . For each channel, 
the receiver was calibrated for differential gain and phase 
offset between the feeds. 

We calculated TOAs from both the WAPP and ASP data 
by using a matc hed-filter algo rithm implemented in the 
Fourier domain ()TavloHri992f) . In this approach, TOAs 
are estimated by comparing the observed total intensity 
pulse profile I {(f)) to a standard template T((?!)). The 
matched-filtering approach is based on the assumption 
that the measured pulse profile is a scaled version of the 
template added to a measurement error term, 

I{^)^aT{cp-cpo) + h + n{4>), (1) 

where a, b and (po are unknown parameters, and n{(j)) is 
noise. Under the assumptions underlying Equation ([!]), 
the TO A uncertainty is P5(j)o, where 5(j)o is the standard 
error of 00- 

For the ASP observations, templates for each ob- 
serving band were generated from high-S/N obser- 
vations of the pulsar co-added from the long-term 
Arecibo/NANOGra\E| timing program, which at present 
uses backend instrumentation and observing bands iden- 
tical to those used in the work presented here (M. Gon- 
zalez, private communication). For the WAPP observa- 
tions, we produced an analytic template containing five 
Gaussian components, consistent with the average pro- 
files in these observations, and c omparable to a previo us 
model of the profile described in lKramer et al.l (|l998f ). 

Residual TOAs were then calculated by fitting a tim- 
ing model ephemeris to the TOAs. For the observations 
here, we used an initial timing ephemeris derived from 
long term monitorin g of J1713-I-0747 cond ucted at the 
Arecibo observatory (jPemorest et al.ll2012[ ). 

3. EVIDENCE FOR PULSE PHASE JITTER 

3.1. Single Pulse Profile Variations 

In this section, we examine single pulses from 
PSR J1713-I-0747. The 50fc individual pulses obtained 
over four minutes allow pulse shape variations to be 
characterized during an interval when the modulation 
from diffractive interstellar scintillation (DISS) is ap- 
proximately constant, since the characteristic DISS time 
scale is about one hour at 1.5 GHz. 

The single pulse signal to noise ratio {S/N, defined as 
the ratio of peak intensity to the rms off-pulse intensity) 
ranged from approximately 1 to 30. In Figure [U we show 
four relatively bright single pulses that have different 
morphologies with peak locations and pulse widths vary- 
ing by amounts that cannot be attributed to radiometer 
noise. 

Figure [2] shows average profiles that have been calculated 
using single pulses selected from different S/N ranges; 
each profile has been normalized to reflect the total de- 
tected flux in the S/N range. Despite occurring three 
times less frequently than the faintest pulses, the bright- 
est and narrowest pulses contribute a total flux that is 

^ www.naiiograv.org 
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Figure 1. : Pulse intensity versus phase (measured in 
cycles, cy) for four single pulses from PSR J1713-I-0747. 
The pulses have been scaled to the same peak amplitude 
to show differences in shape and centroid pulse phase. 

comparable to that of the faintest pulses. We also flnd 
that brighter pulses are narrower with peaks concen- 
trated closer to the leading edge of the average pulse, 
suggesting that brighter pulses have earlier arrival times, 
which we will show below. 

The presence of bright single pulses on the leading edge 
of the integrated pulse can be verified through a num- 
ber of other diagnostics. One of these is the intensity 
modulation index, defined as 



where I{(t>) and (Tj((j)) are the mean and rms intensity 
at phase (j> and cr/. off is the rms off pulse intensity. The 
modulation index measures the normalized rms intensity 
across pulse phase (f) and is sensitive to the frequency of 
occurrence of pulses at each pulse phase. A comparison 
of TO/ {(j)) with the integrated profile (Figure [3]) indicates 
that the modulation index is largest at the leading edge 
of the pulse. This suggests the existence of a tail of 
bright pulses on the leading edge that has large enough 
amplitudes to dominate the modulation index. 

The separations between strong pulses in units of the 
pulse period appear to be random. In Figure |4l we show 
the histogram for the spacing in arrival time AP for 
bright pulses with S/N > 14. The histogram of AP 
is fit well by an exponential distribution and is there- 
fore consistent with a Poisson process. The mean spac- 
ing is ~ 150 pulses. The exponential distribution is the 
same as found for gi ant pulses from the Crab pulsar 
(jLundgren et al.lll995[) . 
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Figure 2. : Average pulse profiles for different ranges of 
S/N. The pulse profiles were formed from single pulses 
with S/N in the different ranges. The ranges are S/N < 
4, 4 < S/N < 10 (labeled on the plot as 4, 10), S/N > 10, 
and for all S/N ranges, and are labeled on the plot at 
a height close to the peak flux for each pulse profile, in 
a line thickness corresponding the thickness of the pulse 
profile. The profiles are normalized to the total fiuence 
so the profiles reflect the relative flux being emitted in 
these S/N ranges. 
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Figure 3. : Top panel: Modulation index to/ as a func- 
tion of pulse phase <j> calculated using equation ([2]). The 
leading edge of the pulse shows the largest rrij. Bottom 
panel: Average pulse profile. 
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Figure 4. : Histogram of the spacing AP of bright 
pulses. The thin line shows the expected exponential 
distribution for AP if the number of bright pulses per 
unit time is Poisson distributed. 



3.2. Single Pulse Timing 

There is a modest correlation between single pulse S/N 
and pulse arrival time. In Figure [5l panel c, we show a 
scatter plot of residual arrival times At versus S/N for all 
pulses with S/N > 3. Brighter pulses tend to have earlier 
arrival times. In panels a and b, we display histograms 
of At for bright {S/N > 20) and weak (14 < S/N < 20) 
pulses, respectively. The weaker pulses arrive in a wide 
range of At, whereas the brightest pulses arrive in a nar- 
rower range of arrival time. While there is a clear cor- 
relation between arrival time and S /N, the correlation is 
not strong and the scatter in arrival times (labeled crj on 
the plot) is comparable in magnitude of the correlation. 
We note that a correlation between S/N and arrival time 
had previously been observed in the bright canoni cal pul- 
sar PSR B0329-h54 (jMcKinnon fc HankinsI [l99l . This 
correlation was identified with mode changing and has a 
different character than the fast jitter effect observed in 
PSR J1713+0747. 

In Figure [5l panel d, we show the histogram for S/N, 
which is a proxy for pulse intensity. For S/N > 3.5, the 
distribution is well described by a log-normal distribu- 
tion. We fit a log-normal distribution to the histogram 
for S/N >3.5 and found a reduced of 1.7. This distri- 
bution is consistent with what is observed in canonical 
pulsars (Cairns et al. 2004). The distribution diverges 
from a log-normal distribution when S/N < 3 because 
radiometer noise (and not intrinsic brightness) dominate 
the estimate of pulse intensity. The distribution is highly 
inconsistent with a power-law distribution because the 
observed distribution does not show a long tail of pulses 
at high S/N. The best-fit power law produced poorly 
matched the data and had a reduced of 150. 



3.3. Timing Fluctuations for Integrated Profiles 
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Figure 5. : Single pulse amplitude and arrival time 
statistics for PSR J1713+0747. Panel a: Histogram 
of the arrival times for very bright {S/N > 20) sin- 
gle pulses. Panel b: Histogram for weaker single pulses 
(14 < S/N < 20). Panel c: Scatter plot of residual ar- 
rival time At versus S/N for single pulses. The horizon- 
tal dashed lines indicate the S/N ranges used for panels 
a and b. The vertical dashed-dot line shows the mean 
pulse arrival time. For clarity, only a random selection 
of the single pulses have been plotted. To compensate for 
fewer bright single pulses, the probability of a pulse be- 
ing selected was scaled oc 1 / p{S / Ni) , where p{S/N) is the 
probability density function for S/N derived empirically 
from panel d. Panel d: Histogram of the S/N of single 
pulses. On the histogram, we also plot the Icr counting 
errors. The best-fit curve to a log-normal distribution is 
displayed as a thin line. 



Analysis of integrated pulse profiles produced by the ASP 
instrument confirms that intrinsic pulse shape variations 
are manifested as TOA variations in the integrated pro- 
files. 

In Figure [HI the residual TO As for ASP observations are 
displayed for observations from MJD 54978, a day in 
which diffractive and refractive scintillation enhanced the 
flux in the 64 MHz ASP band by a factor of « 4. Within 
the figure, we show residual arrival times for simultane- 
ous observations within three sub-bands of 4 MHz band- 
width, widely spaced across the band. The error bars 
represent the white noise uncertainty in the TOA esti- 
mates and vary in size between the sub-bands and in time 
because the flux (though relatively high in all bands) is 
modulated by diffractive scintillation. 

Figure [6] also shows histograms of the residual TOAs for 
the three sub-bands. The values for skewness and kur- 
tosis that we estimate indicate that the histograms are 
consistent with Gaussian distributions. 

The arrival times are strongly correlated between the 
sub-bands. In Figure [3 we show a scatter plot of 
the residual arrival times in two frequency bands. A 



best fit line to the data shows a slope consistent with 
unity, implying that the white-noise component of the 
TOA variations not caused by radiometer noise is es- 
sentially independent of frequency. This analysis was 
conducted for all pairs of frequency channels, and in all 
cases a slope consistent with unity was found. We note 
that a similar correlation between sub-banded TOAs was 
found in Park e s obs ervations of PSR J0437— 4715 by 
lOslowski et al.l ()2011f ). suggesting that the effect is not 
observatory dependent. 
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Figure 6. : Left panels: Residual TOAs for obser- 
vations simultaneously obtained in three sub-bands of 
4 MHz bandwidth. The central frequencies of the bands 
are 1382 MHz for panel a, 1410 MHz for panel 6, and 
1430 MHz for panel c. Right panel: Histograms of resid- 
ual TOAs. 



Another demonstration of the large correlation of the 
TOA variations across frequency is given in Figure |S1 
which shows the cross correlation functions (CCFs) for 
arrival time variations between the three pairs of sub- 
bands. The CCFs show spikes at zero lag that drop to 
zero correlation after one sample (10 s). The zero-lag 
spike confirms that the residual TOAs are correlated be- 
tween different frequency channels, and that the struc- 
ture is unresolved in time at the 10 s level. This is con- 
sistent with what would be expected if the single pulse 
modulation affecting the WAPP observations was caus- 
ing the timing error observed in the ASP observations, 
because the variations appear to be uncorrelated with 
time. 

If we define the rms pulse-phase jitter aj as the excess 
white-noise residual over that expected from radiome- 
ter noise, we can estimate it from the CCF because ra- 
diometer noise will not correlate between different fre- 
quency bands while the jitter will. The zero-lag value 
is CCF(O) = crj ~ 0.30 Ais^, so crj w 0.55 pa for T = 10 s 
(2188 pulses). This corresponds to jitter in single-pulse 
arrival times of -\/2188 x 0.55 ps w 26 ps, compared to 
the pulse width (FWHM) of 110 ps. Our observation of 
shape variations are consistent with a model presented 
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Figure 7. : Correlation of TOAs between two frequency 
channels. The error bars represent the \ — a template 
fitting errors. The dashed angled line indicates ^ty = 

in iCordes fc Downj ()1985[ ). in which single pulses are 
modeled to have narrower profiles than the mean pro- 
file, but large variations in centroid. Shape variations 
can then be characterized by a dimensionless jitter pa- 
rameter fj — 1 — (wi/wa)'^, where Wi is the width of a 
single pulse and Wa is the width of an average pulse. For 
PSR J1713-^0747, we find a jitter parameter fj « 0.24, 
which is consistent wit h values found for other pulsars 
(|Cordes fc DownsllT985l) . 

3.4. Connecting Single Pulse Timing Variations to 
Standard Precision Timing Observations 

Phenomenologically, it has been known since the early 
days of pulsar timing that despite the phase jitter 
of individual pulses, integrated profiles converge to- 
ward an average shape that is stable on time scales of 
decades. This convergence corresponds to timing pre- 
cision that improves as 1/ ^/N^ when Np pulses are 
summed (jHelfand et aLlflOTSl ) . For some long-period pul- 
sars, the effective number of independent pulses is less 
than Np because there are short term correlations and 
the results for J1713+0747 are consistent with a lack of 
correlations. 

We demonstrate that the arrival-time errors from stan- 
dard template fitting to average profiles obtained by av- 
eraging Np ^ 1 pulses are completely consistent with the 
directly measured pulse-to-pulse variations. We do so by 
comparing the actual timing errors from ASP data with 
the errors expected under idealized conditions where an 
average profile is the sum of the template and additive 
white noise. In the top panel of Figure [HI the rms resid- 
ual o'iyNp) is plotted against the number of pulses aver- 
aged, Np, for both ASP and WAPP observations (open 
symbols). We also show rms timing residuals ar{Np) ex- 
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Figure 8. : Cross correlation function between time se- 
ries for the 1430 MHz channel and three other channels. 
The lower dashed lines show the zero correlation level. 
The CCF shows a spike at zero lag that decays to zero 
in one sample (10 s). The upper dashed lines shows the 
level of jitter inferred from Figure |9l 

pected if there were only additive radiometer noise and 
no phase jitter (filled symbols). Values for cr^ were esti- 
mated from simulated pulse profiles calculated by adding 
white noise to the template shape using the S/N appro- 
priate for the data. The simulated profiles were ana- 
lyzed identically to the observations to produce arrival 
time estimates. These simulated profiles also satisfy the 
conditions for matched filtering and therefore yield the 
smallest possible timing error under the assumption that 
there is minimally staMonary (white) noise present in the 
timing observations (Turin 1960). The predicted rms er- 
ror for additive noise alone is much less th an a. Both a 
and CTr decrease proportionally to l/^^Np, as expected 
for timing errors that are uncorrelated between pulses. 

In the bottom panel of Figure [9l we show the quadrature 
difference between the actual and predicted timing errors 

aj{Np) = ^aHNp)-a^ANp). (3) 

After accounting for the different numbers of pulses av- 
eraged, the values of aj obtained from the WAPP and 
ASP observations consistent with each other. The single- 
pulse value 'Jj(l) ~ 27± 1 /is is consistent with the value 
derived from the cross-correlation analysis in the previ- 
ous section. The total timing error is a combination of 
both phase jitter and radiometer noise. 

4. CONFIRMING THE JITTER HYPOTHESIS 

We have shown that the excess timing error is broadband 
and uncorrelated between pulse profiles calculated from 
disjoint sets of pulses. In addition to jitter, there are two 
other plausible forms of timing error that can contribute 
errors of this nature: diffractive interstellar scintillation 
and observation calibration. Here we rule out both of 
these alternative explanations. 
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Figure 9. : Top panel: Measured rms (unfilled symbols) 
and expected rms (filled symbols) timing residuals versus 
the number of pulses averaged Np, for observations using 
the WAPP (triangles) and ASP backends (circles). Bot- 
tom panel: Quadrature difference between observed and 
expected rms timing errors. The dashed line is a fit to 

— 1/2 

the data, showing the expected scaling (Tj(iVp) c>c iVp 
Despite having different levels of radiometer noise, the 
level of jitter extrapolated from the WAPP observations 
is consistent with the ASP observations. The horizontal 
scale at the top of the plot gives the integration time in 
seconds of sub-integrations Aisub corresponding to Np in 
the bottom scale. 



Diffractive interstellar scintillation (DISS) will generally 
modulate the pulse shape on minutes to hour time scales, 
depending on observing frequency, dispersion measure, 
and direction. The shape perturbation is correlated 
over a time equal to the diffractive time scale td and 
over a frequency range equal to the scintillat ion band- 
widthj t/rf (Cordes ct al. 1990; Hemberger fc Stinebring 
[200a ICordes fc ShannonI I20iai The DISS time scale 
and bandwidth vary strongly with observing frequency, 
approximately as and i'"*, respectively. The DISS 
timing perturbation has been observed in another (but 
much more strong l y scattered) MSP , PSR BI937-I-21 
(|Cordes et al.l (19901: iJenet fc Gi]|l2004D . 

DISS is ruled out as the source of the excess TOA vari- 
ation observed for MSP J1713-K0747. First, diffractive 
TOA variations would be correlated over a time compa- 
rable to the 1 hr diffractive time scale whereas the TOA 
variations appear to vary on a time scale shorter than the 
rotation period of the pulsar (4.6 ms). Also, the timing 
perturbation expected from DISS is much smaller than 
the observed timing error. Diffractive scintillation is ex- 
pected to contribute an rms error of (Cordes ct al...l990i : 



At 



DISS 



(4) 



where r is the pulse broadening time and is the num- 
ber of "scintles," which are the number of bright patches 
of constructive interference contained in the frequency- 
time plane of an observation. The number of scin tles is 
Ns w (l-h0.3T/id)(l + 0.3Ai//Ai/d) (|Cordes et al.l [T990'). 
where T is the observing time for each sub-integration 
(which for our observations vary in duration between the 
rotation period of a pulsar 4.6 ms and « 10^ s), td is 
the diffractive time scale (> 1 hr), Av is the observing 
bandwidth 100 MHz and Ai/d is the diffractive scintilla- 
tion bandwidth (sa 50 MHz). Thus for the observations 
reported here, ~ 1. The pulse broadening time can be 
inferred from the scintillation bandwidth r = Ci/2TTAi>di 
where Avd is the sci ntillation bandwidth an d Ci is a con- 
stant of order unity ([Cordeii&UicketllUll). Evaluating 
Equation ([J), we find that Aioiss ~ 3 ns, which is much 
less than the rms scatter of 26 /is. 

It is also possible that instrumental effects (including 
both calibration errors and radio frequency interference) 
can cause pulse shape changes. However we disfavor 
this origin. We convincingly connect the levels of jit- 
ter observed in ASP observations with that observed in 
WAPP observations. These backends operate with dif- 
ferent architectures. They also have different sampling 
levels, suggesting that artifacts of digitization , such as 
scattered power (jKouwenhoven fc Voutel[200lf ). are not 
grossly affecting the 3-level sampled WAPP observations. 
The jitter is observed at modestly different frequencies 
(ruling out narrow-band RFI). In a previous observing 
campaign in which PSR J1713-I-0747 was observed si- 
multaneously with the Arecibo and Green Bank tele- 
scopes, residual TOAs were found to be correlated on 
similarly short time scales at t he same observing fre- 
quency (|Lam fc Demores"^ 120101 ). Calibration errors as- 
sociated with non-orthogonality of the feeds should be 
uncorrelated between the different telescope feeds. Ad- 
ditionally, timing errors induced will change slowly as the 
orientation of the feed with respect to the pulsar changes. 
Radio frequency interference (RFI) could also plausibly 
cause timing error. Interference is local to each telescope 
and will be uncorrelated between telescopes separated by 
thousands of kilometers. 

The correlation between larger S/N and earlier arrival 
times is seen only at the single-pulse level and is therefore 
only observed in the WAPP data sets. Our interpretation 
is that the effect is evidently diminished in averages of 
10 s or longer. There is no evidence that the correlation 
is caused by instrumental effects. Firstly, the effects of 
digitization are small because the signal to noise ratio of 
the detected, non-dispersed signal is low. The dispersion 
smearing across the band is approximately 3 ms (much 
greater than the pulse width) so the non-dispersed signal 
has an S/N of approximately 0.7. Distortions in pulse 
shapes are strongest when the dispersive delay across 
the band is comparable to the pulse width and are sup- 
pressed wlien^_tlie delay is much larger than the pulse 
width (jStairs et al.ll2000l) . The dominant effect of dig- 
itization is the appearance of negative detected power 
(dips) at the edges of the pulse profile (Jenet fc Andersolil 
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Shannon & Cordes 



[1991 . We see no evidence for this type of distortion in 
our data. For example in Figure [2l the basehne is flat 
at the edge of even the brightest pulses {S/N > 10) that 
would be the most affected by distortions. The most con- 
vincing demonstration of an intrinsic origin for the corre- 
lation between S/N and arrival time would be to observe 
the same effect with different instrumentation that pro- 
vides fast-sampled data having more bits of representa- 
tion (e.g., > 8 bits) and preferably different architecture 
(e.g., a baseband recorder or a polyphase filterbank). 

5. CORRECTING TOAS FOR PULSE SHAPE CHANGES 

We investigated two distinct methods for correcting 
arrival times due to short term variations: multi- 
component template fitting and principal component 
analysis (PCA). These techniques were tested on both 
ASP and WAPP data, enabling us to test the correction 
algorithms on data over a wide range of pulse averaging. 
We were unable to develop an algorithm to improve the 
precision of arrival times using pulse shape information. 
We attribute the lack of success to the absence of a suf- 
ficiently strong correlation between pulse shape changes 
and arrival times presented in the previous sections. Here 
we briefly summarize two attempts at correction and de- 
fer further discussion to a future paper. 

Wc implemented a multi-component TOA fitting algo- 
rithm similar to one that was u sed to improve the timing 
precision of PSR J1022-f-1001 ()Kramer et al.ill999il . For 
PSR J1713-f0747, we constructed an analytic template 
containing 5 Gaussian co mponents, consistent with pre- 
vious models of the pulsar ([Kramer ei~aII[T998h . We then 
made a non- linear fit for the shape parameters, i.e., the 
amplitude, widths, and relative positions of the different 
components. The TOAs generated with this technique 
were then compared to TOAs inferred using the stan- 
dard method. We found no improvement in the resid- 
ual arrival times from TOAs derived from the non-linear 
technique. 

We also conducted principal component analysis (PCA) 
on our data using a meth od comparable to one presented 
in lOslowski et al.l ([201 1[ ) , though implemented using a 
different algorithm. 

To bricfiy summarize our algorithm, we calculated data 
vectors by subtracting the average pulse profile from in- 
dividual profiles. The covariance matrix of the data vec- 
tors was then diagonalizcd to obtain the eigenvalues and 
eigenvectors. 

A scalar quantity, formed by projecting the observed 
pulse profiles onto the most significant eigenvectors, is 
then calculated and correlated with the residual TOAs. 
Any significant correlation between the scalar and resid- 
ual TOAs indicates that a particular eigenvector con- 
tains aspects of profile variati on that is influencing arrival 
times. Oslow ski et al.l ([201 ID found such an eigenvector 
for PSR J0437— 4715 and were able to use the correlation 
to improve the timing precision of the pulsar. For PSR 
J1713-f 0747, we found no significant correlation between 
any eigenvector and residual arrival time and were there- 
fore unable to implement any correction technique using 
TOAs. 

For both PSR J1022+1001 and PSR J0437-4715, cor- 



related pulse shape changes are highly significant in in- 
tegrated pulse profiles, with sub-integrations separated 
by minutes to days showing differences in pulse shape 
visible by eye. For PSR J1713-I-0747, the pulse shape 
changes appear to decorrelate on pulse to pulse time 
scales on the millisecond-hour baselines probed by the 
observations presented here. 

6. DISCUSSION AND CONCLUSIONS 

We have demonstrated that pulse jitter is the dom- 
inant source of timing error in the observations of 
PSR J1713-I-0747 presented here. This pulsar is cur- 
rently being monitored by the three major pulsar tim- 
ing array collaborations: the European Pulsar Timing 
Array (EPTA, Fcrdman ct al. 2010), the North Amer- 
ican Nanoher tz Observatory fo r Gravitational Waves 
(NANOGrav, iJenet et al.[|2009D. and the Pa rkes Pulsar 
Timing Array (PPTA. iVerbiest et al.|[2Q10l ). Based on 
the observations presented here, the timing error induced 
by shape variations can be extrapolated to the longer ob- 
serving spans typically used in pulsar timing array ob- 
servations. For a 30 min observation with Np w 10^'^, 
we predict aj « 40 ns. This is consistent with the rms 
error o bserved in the long ter m timing precision of this 
object (iDemorest et al.l[2012t ). suggesting that a signif- 
icant component of the timing error is associated with 
jitter. 

Pulse shape changes, if uncorrelated on time scales longer 
than typical observing cadences and uncorrectable using 
information about pulse shape changes, contribute addi- 
tional white noise to TOAs. In general there are three av- 
enues for reducing white noise in timing observations: ob- 
serving with a more sensitive telescopes, observing with 
a wider observing bandwidth, or averaging over a larger 
number of pulses. Of these three options, only the last 
can be used to reduce timing error associated with jit- 
ter because the TOA perturbations are broadband and 
independent of antenna gain. This implies that many 
pulsars may need to be observed with higher throughput 
(i.e, with longer observations at each epoch, or observed 
with shorter cadence) to make a convincing detect ion of 
the g ravitational wave background (Cordes fc Shannon! 
[2011 . 

The timing precision of PSR J1713+0747 and probably 
other bright MSPs is likely limited by jitter if observed 
with telescopes where the single pulse S/N > 1. If we 
assume that the pulse is Gaussia n shaped with width w 
the rms associated with j itter is ([Cordes fc Downslll985l : 
ICordes fc Shannon! [2011 

Atj^wJj{l + mj)N~'/^, (5) 

where Wi is the intrinsic pulse width, fj is the jitter 
parameter, m/ is the modulation index, and N is the 
number of pulses used to form each subintegration. In 
contrast, the rms expected from radiometer noise is 



where Wa is the width of the integrated pulse and 5*0 
is the single pulse S/N. We can solve for the critical 5*0 
by equating equations (O and ([6]). We find that the 
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threshold single pulse signal to noise ratio is 



(7) 



/,/(!- /,7)i/2(l + m2) 



1.4 



0.5 



0.5 



The exact value of the transition depends on the magni- 
tude of phase jitter, the modulation index, and the shape 
of the pulse. However we expect that for most pulsars 
the critical single pulse S/N will be close to unity. 

Large telescopes provide other avenues for improving 
TOA precision. First, higher S/N observations will be 
able to better diagnose the pulse shape changes in a 
larger number of MSPs. Additionally, larger interfero- 
metric telescopes such as the SKA can provide a higher 
timing throughput if they possess the ability to form mul- 
tiple sub-arrays. In this case, multiple pulsars can be 
observed simultaneously, with the gain of each sub-array 
tailored to individual pulsars. 

Timing errors associated with shape changes are cor- 
rectable only if they are accompanied by gross changes in 
the shape of the pulse. If not, correction is likely minimal 
because other astrophysical effects will cause similar dis- 
placement of the pulses (albeit on longer time scales). 
Though unsuccessful in the analysis of J1713-I-0747, 
multi-component template fitting and PCA are sensitive 
to subtle changes in the shape of pulses. 

It is imperative to assess profile stability in all current 
and candidate MSPs for pulsar timing arrays. As we 
have demonstrated, single pulse measurements are not 
required to determine the amount of pulse phase jitter. 
It can readily be inferred from integrated profiles. 
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